On convergence to equilibrium distribution 

for Dirac equation 



A. Komech0' 2 

Faculty of Mathematics of Vienna University, 1090 Vienna, Austria 
and IITP RAS, Moscow, B.Karetny, 19 
e-mail: alexander.komech@mat.univie.ac.at 

E. Kopylovad 
IITP RAS, Moscow, B.Karetny, 19 
e-mail : elena . kopy lova@univie . ac . at 



Abstract 

We consider the Dirac equation in R 3 with a potential, and study the distribution \i t 
of the random solution at time t € R. The initial measure ^jlq has zero mean, a translation- 
invariant covariance, and a finite mean charge density. We also assume that uq satisfies 
a mixing condition of Rosenblatt- or Ibragimov-Linnik-type. The main result is the long 
time convergence of projection of m onto the continuous spectral space. The limiting 
measure is Gaussian. 
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1 Introduction 



This paper can be considered as a continuation of our papers [5]- [Hj, [H] which concern the long 
time convergence to equilibrium distribution for the linear wave, Klein-Gordon and Schrodinger 
equations. 

The convergence should clarify the distinguished role of the canonical Maxwell-Boltzmann- 
Gibbs equilibrium distribution in statistical physics. One of fundamental examples is the 
Kirchhoff-Planck black body radiation law which specify the equilibrium distribution for the 
Maxwell equations, and served as a basis for creation of quantum mechanics. The law likely 
should be correlation function of limiting equilibrium measure for coupled Maxwell- Schrodinger 
or Maxwell-Dirac equations. 

Our ultimate goal would be the proof of the convergence for nonlinear hyperbolic PDEs. At 
the moment, a unique result in nonlinear case has been proved by Jaksic and Pillet for wave 
equation coupled to a nonlinear finite dimensional Hamiltonian system |12j . 

The main peculiarity of the problem is the time-reversibility of dynamical equations. For 
infinite particle systems this difficulty was discussed in Boltzmann-Zermelo debates (1896-1897). 
Many attempts were made to deduce the convergence from an ergodicity for such systems by 
H. Poincare, G. Birkhoff, A. Hinchin, and many others. However, the ergodicity is not proved 
until now. 

In 1980 R. Dobrushin and Yu. Suhov have introduced a totally new idea for obtaining the 
convergence to equilibrium measures imposing a mixing condition on initial distributions [I] in 
the context of infinite particle systems. 

We develop this approach for hyperbolic PDEs. In [5]- [8], [10]- [11] the convergence to 
equilibrium distributions has been proved for the linear wave, Klein-Gordon and Schrodinger 
equations with potentials, for the harmonic crystal, and for the free Dirac equation. The 
initial distribution are translation invariant and satisfy the mixing condition of Rosenblatt or 
Ibragimov-Linnik type. 

Here we consider the linear Dirac equation with the Maxwell potentials in M. 3 : 



where ip(x, t) e C 4 , m > and a = (ai, a 2 , a 3 ). The hermitian matrices (3 = a and satisfy 
the following relations: 



iip(x, t) = Hip{x, t) 
4>(x,0) = ip Q (x) 



[-ia ■ V + f3m + V(x)]if)(x, t) 



xem 3 



(l.i) 




k, I = 0,1,2,3,4. 



The standard form of the Dirac matrices and j3 (in 2 x 2 blocks) is 




(1.2) 



where I denotes the unit matrix, and 




(1.3) 
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We assume the following conditions: 

El. The potential V G C°°(1R 3 ) is a hermitian 4x4 matrix function such that 

\d a V(x)\ < C{a)(x)- p - lal , (x) a = (1 + \x\ 2 ) a/2 (1.4) 

with some p > 5. 

E2. The operator H presents neither resonance nor eigenvalue at thresholds. 

Under the condition E2 the operator H has a finite set of eigenvalues Uj G (— m, m), j = 1, N 

1 k ■ 

with the corresponding eigenfunctions Q,...,Q J , where kj is the multiplicity of ujj. Denote by 
Pj the Riesz projection onto the corresponding eigenspaces and by 

P c :=l-P d , P d = Y< P i ( L5 ) 

3 

the projections onto the continuous and discrete spectral spaces of H. 

We fix an arbitrary 5 > such that 5 + 5 < p and consider the solutions ip(x, t) G C 4 with 
initial data i[>o( x ) which are supposed to be a random element of the weighted Sobolev space 
"H = L 2 _ 5 / 2 _ s , see Definition 12.11 below. The distribution of ip is a Borel probability measure 
p on % with zero mean satisfying some additional assumptions, see Conditions S1-S3 below. 
Denote by pt, t G R, the measure on 7i, giving the distribution of the random solution ip(t) to 
problem ( II. ip . We identify the complex and real spaces C 4 = M 8 , and <g> stands for the tensor 
product of real vectors. The correlation functions of the initial measure are supposed to be 
translation-invariant: 

Qo(x,y) :=E(i/) (x)®i/) (y)) =q {x-y), x,y E 1R 3 . (1.6) 

We also assume that the initial mean charge density is finite: 

e := E\ij (x)\ 2 = trg (0) < oo, x G 1R 3 . (1.7) 

Finally, we assume that the measure p$ satisfies a mixing condition of a Rosenblatt- or Ibragimov- 
Linnik type, which means that 

ipo(x) and ipo{y) are asymptotically independent as \x — y\ — > oo. (1-8) 

Let P*pt denote the projection of p t onto the space H c := P c %. Our main result is the (weak) 
convergence of P*pt to a limiting measure u^, 

P*p t —r Voo, t ->■ oo, (1.9) 

which is an equilibrium Gaussian measure on H c . A similar convergence holds for t — > — oo 
since our system is time-reversible. 

The convergence ( II .9p for the free Dirac equation with V(x) = has been proved in [8]. 
The case of the perturbed Dirac equation with 7^0 requires new constructions due to the 
absence an explicit formula for the solution. To reduce the case of perturbed equation to the 
case of free equation we formally need a scattering theory for the solutions of infinite global 
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charge. We manage a dual scattering theory for finite charge solutions to avoid the infinite 
charge scattering theory: 

P c U'{t)4> = U' Q {t)W4> + r(t)<j>, t > 0. (1.10) 

Here U' Q (t) and U'(t) are a 'formal adjoint' to the dynamical groups Uo(t) and U(t) of the free 
equation with V = and equation ( 11. ip with V ^ respectively. The remainder r(t) is small 
in the mean: 

£|(^ o ,r(t)0)| 2 ^O, t-»oo. (1.11) 

where (•, •) is defined in (I2.20p . This version of scattering theory is based on the weighted 
energy decay established in [2]. 

2 Main results 
2.1 Well posedness 

Definition 2.1. For s, a 6 R, let us denote by H% = i/*(lR 3 ,C 4 ) the weighted Sobolev spaces 
with the finite norms 

\ H s = \\(xr(vm\ L2 <oo. 



We set L? a = H®. Note, that the multiplication by V(x) is bounded operator L 2 — > L 2 +p . 
The finite speed of propagation for equation fll.ip implies 

Proposition 2.2. i) For any ip G L?_ a with < o < p there exists a unique solution ip(-,t) G 
C(R, L 2 __ a ) to the Cauchy problem (fl.lj) . 

ii) For any i el, the operator U(t) : ip Q \-t ip(-,t) is continuous in L 2 _ a . 



Proof. Fist, consider the free Dirac equation: 



x(x, t) = Hqx{ x -> t) = {—ot ■ V — if3m)x(x, t) x G M 3 , 
x(x,o) = i>o(x). 



(2.12) 



Let s6R and ipo G L 2 . In the Fourier space the solution to ( 12. 12ft reads 

x (k,t) = e i(a - k - pm)t ip (k). 
Since ipo £ H s then x{-,t) G and the bounds hold 

llx(-,*)lk = *)!!*• < < ciWIIV^Ik • ( 2 - 13 ) 



Now consider perturbed equation (II. ip . Let < a < p and ip G L 2 _ a . We seek the solution to 
(11.11) in the form 

i/)(x,t) = X (x,t) + (/>(x,t), (2.14) 
where x(t) = U (t)ip G L?_ a is the solution to free equation (12.121) . and 

<p(x,t) = H(p(x,t) + Vx(x,t), 0(x,O) = O. (2.15) 
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Since 0(0) = and V\ G L 2 then there exists the unique solution 0(t) G L 2 to (12.151) which is 
given by Duhamel representation: 

t 

0(t) = ju{t- r)V X (r)dr. 
o 

Finally, by charge conservation for the Dirac equation we obtain 

\\U(t-r)V X (r)\y a < \\U{t-r)V X {r)\y = \\V X (t)\\& < C\\ X {r)\\ L% < C7|| X (r)||^ < oo . 

□ 



2.2 Random solution. Convergence to equilibrium 

Let (Q, E, P) be a probability space with expectation E and B(H) denote the Borel cr-algebra 
in "H. We assume that 0>o = ipo(u),-) in (11.11) is a measurable random function with values 
in ("H, B(H)). In other words, (u,x) H- i[)q(cu,x) is a measurable map Q x ]R 3 — y C 4 with 
respect to the (completed) cr-algebras £ x i3(IR 3 ) and i3(C 4 ). Then, owing to Proposition 12.21 
%fj{t) = U{t)ipo is again a measurable random function with values in (%,&(%)). We denote 
by no{dipo) a Borel probability measure in 7-L giving the distribution of the random function 
ipQ. Without loss of generality, we assume (f2,£,P) = B(T-L), fi,) and ipo(u,x) = oj(x) for 
fi (du) x rfx-almost all (u,x) G H x M. 3 . 

Definition 2.3. fit is a probability measure on "H which gives the distribution of ip{t): 

fi t (B) = iM (U(-t)B), V.B G B(H), t>0. (2.16) 

Denote by P*^t the projection of measure /i f onto % c = P c l-L: 

P: f i t (B) = fi t (P c - 1 B), VBeB(U c ), t>0. (2.17) 

Our main goal is to derive the weak convergence of P*[J>t in the Hilbert space P c HZl for any 
e > 0, and a > 5/2 + 5: 

P*Ht — ? as t — > oo, (2-18) 

where is a Borel probability measure on P C HZ^. By definition, this means the convergence 

J /(^)P c */x*(#) -> y /(V>)foo(#) as t -> oo. (2.19) 

for any bounded and continuous functional /(VO i n Pc.B-%- 

Set 72/0 = (R- e ^5 Im-0) = {Re • • • , Re '04, Im0i, . . . , Im-0 4 } for ip = (ipx, . . . -0 4 ) G C 4 and 
denote by IV ip the j-th component of the vector TZip, j = 1, 8. The brackets (•, •) mean the 
inner product in the real Hilbert spaces L 2 = L 2 (R 3 ), in L 2 (g) ~R N , or in some their different 
extensions. For ip(x),<f)(x) G L 2 (M 3 ,C 4 ), write 

8 

(0,0) := (7^,^0) = ]T (7^0, 7^0). (2.20) 
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Definition 2.4. The correlation functions of the measure fi® are defined by 

Q%(x,y) = S^Vo(aj)^- J Vo(y)) for almost all x,y E M 3 , i,j = 1,...,8, (2.21) 

provided that the expectations in the right-hand side are finite. 

Denote by D the space of complex- valued functions in C(f(M 3 ) and write T> := [D] 4 . For a 
Borel probability measure fx denote by fx the characteristic functional (the Fourier transform) 

P>{<j>) = J exp(z(^,0))/i(#), G V. 

A measure /x is said to be Gaussian (with zero expectation) if its characteristic functional is of 
the form 

fi{4>) =exp{-^Q(0,0)}, (ft eT>, 

where Q is a real nonnegative quadratic form onD. A measure /i on "H is said to be translation- 
invariant if 

Li(T h B) = Li(B), BeB(H), h G R 3 , 
where T h ijj(x) = if>(x — h), x G M 3 . 

2.3 Mixing condition 

Let O(r) be the set of all pairs of open bounded subsets A, B C 1R 3 at the distance not less 
than r, dist(*4, B) > r, and let a (A) be the a-algebra in H generated by the linear functionals 
ip i — y (V'>0)> where <p E T> with supp0 C A. Define the Ibragimov-Linnik mixing coefficient of 
a probability measure llq on % by the rule (cf. [9] Def. 17.2.2]) 

/ n _ \Li Q (AnB)- li (A)ix (B)\ 

(p[r) = sup sup — . (2.22) 

(A,B)&o(r) a e a (A), B G a(B) ^ B > 

Vo(B) > 

Definition 2.5. We say that the measure llq satisfies the strong uniform Ibragimov-Linnik 
mixing condition if 

ip(r) ^ as r -»■ oo. (2.23) 
We specify the rate of decay of cp below (see Condition S3). 

2.4 Main assumptions and results 

We assume that the measure fx has the following properties SO— S3: 

SO fx has zero expectation value, 

Eijj Q {x) = 0, x G M 3 . 
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51 fiQ has translation invariant correlation functions, 

Q%(x,y) = E(K i i> Q {x)K j My))=<li j (x-y), £,J = 1, 8 (2.24) 
for almost all x, y G M 3 . 

52 /i has finite mean charge density, i.e. Eqn (jl.7p holds. 

53 /io satisfies the strong uniform Ibragimov-Linnik mixing condition, with 

r 2 ip 1/2 (r)dr < oo. (2.25) 

Remark 2.6. The examples of measures on L 2 oc (R 3 ) satisfying properties S0-S3 have been 
constructed in JE/ (see $$2.6.1-2.6.2). The measures on L 2 _ a with any a > 3/2 can be construct 
similarly. 

Introduce the following 8x8 real valued matrices (in 4 x 4 blocks) 

Note that = A k , k = 1, 2, 3, = -A . Write 

A = (A 1 ,A 2 ,A 3 ), P = A-V + mA . (2.27) 
For almost all x,y £ M 3 , introduce the matrix-valued function 

Qoo(x, y) = (q^(x, y)) = (q%{x - y)) . (2.28) 

V / i,j=l,...,8 V / ij=l,...,8 

Here 

M*) = ^q (k) - ~V(k)P(k)q (k)P(k), (2.29) 

P(/c) = — zA- k + mAo, V{k) = 1/ (k 2 + m 2 ), and go(&) is the Fourier transform of the correlation 
matrix of the measure //q (see 12.24]) . We formally have 

qoo(z) = ~q (z) + l -V * Pq (z)P (2.30) 

where V(z) = e _m ' 2 '/(47r|2;|) is the fundamental solution for the operator —A + m 2 , and * 
stands for the convolution of distributions. 

Lemma 2.7. Let conditions SO, S2 and S3 hold. Then 

q eLP(R 3 ), p>l. (2.31) 

Proof. Conditions SO, S2 and S3 imply (cf. [HI Lemma 17.2.3]) that 

\q^z)\<Ceo<py 2 (\z\), z E R 3 , z,j = l,...,8. 

The mixing coefficient if is bounded, hence 

/roc 
p p/2 (\z\)dz <dj r 2 cp 1/2 (r)dr < oo 

by flZZgp . □ 
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Lemma [2TT1 with p = 2 imply that % E L 2 . Hence, qoo E L 2 by (I2.29p . and q^ also belongs 
to L 2 by (I2T30|) . 

Denote by a real quadratic form on L 2 defined by 

8 

Qoo(0,0) = (QoofayWix)®-*^)) = V / Q^O, y)n l (t){x)n 3 (t)(y)dxdy 

a ^IR 3 xR 3 

Corollary 2.8. The form is continuous on L 2 because qo(k) and then qoo{k) are bounded 
by Lemma WH\ and formula / [2.29|) . 

Our main result is the following: 

Theorem 2.9. Let m > 0, and let conditions El— E2 ; SO— S3 hold. Then 

i) the convergence in / [2.18)) holds for any e > and o > 5/2 + 5. 

ii) the limiting measure /i^ is a Gaussian equilibrium measure on ~K C . 
Hi) the characteristic functional of is of the form 

= exp{~Q oo (W0,W(j))}, (ft E T>, 
where W :T> — >• L 2 is a linear continuous operator. 

2.5 Remark on various mixing conditions for initial measure 

We use the strong uniform Ibragimov-Linnik mixing condition for the simplicity of our presen- 
tation. The uniform Rosenblatt mixing condition [13] with a higher degree > 2 in the bound 
fll.7p is also sufficient. In this case we assume that there exists an e, e > 0, such that 

sup E\%1)q(x)\ 2+£ < oo. 
Then condition (12.251) requires the following modification: 

/■oo 

/ ra p (r)dr < oo, p = min(e/(2 + e), 1/2), 
Jo 

where a(r) is the Rosenblatt mixing coefficient defined as in (I2.22p . but without the denominator 
Ho(B). The statements of Theorem 12.91 and their proofs remain essentially unchanged. 

3 Free Dirac equation 

Here we consider the free Dirac equation (I2.12p We have 

(d t - a ■ V - i/3m)(d t + a ■ V + ifim) = d 2 - A + m 2 
Then the fundamental solution G(x, t) of the free Dirac operator reads 

G t (x) = (d t - a ■ V - i/3m)S t (x) (3.1) 
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where St(x) is the fundamental solution of the Klein-Gordon operator d 2 — A + m 2 : 

£t(x) = F-lf-^, L0 = L0(k) = y/\k\ 2 + m 2 . (3.2) 
00 

Using the notations (12.261) and ( 12. 27ft . we obtain in real form 

Kx(t) = g t *KiPo, Gt = (d t -P)£ t . (3.3) 

The convolution exists since the distribution £t{x) is supported by the ball \x\ < t. Now we 
derive an explicit formula for the correlation function 

Q t (x, y) = q t {x-y) = E (tZ X {x, t) ® TZ X (y, *)) (3.4) 

Lemma 3.1. (cf. |2l Formula (4-6)]) The correlation function Qt(x,y) reads 

1 _ 1 [l + cos2u)t„ ,, s sin2o;t. 



Q t (x,y) = q t (x - y) = F k \ x _ y q (k) - (q (k)P(k) - P(k)q (k)) 

1 — cos 2ut 



2 v ' 2co 

P(k)q {k)P(k)] (3.5) 



2a; 2 

Proof. Applying the Fourier transform to (13. 3 j) we obtain 

TZ^(k,t) =Gt(kyR$ (k) = (coswt-P(Jfe)^^)4(fe) (3-6) 
By translation invariance condition (I2.24p we have 

E{K4> Q {k) ® 7^ (*0) = ^fe^^o^ " y) = (27r) 3 6(k - k')q (k) 
Then ( 13. 6p implies that 

E(Kx~(k,t)®Kx~(k',t)) = {2nf5{k - k')g t {k)q {k)g* t {k) 

Therefore, 

q t (k) = g t (k)q (k)g:(k) = (cosut-P(kf-^)q (k)( C osut + P(kf-^) 

\ 00 J \ OJ / 

since P*(k) = -P(k). Hence ([33]) follows. □ 
Corollary 3.2. For any z G K 3 £/ie convergence holds 

qt{z) ->■ 900(2), t -)> 00 

where q^z) is defined in ( 12.30)) . 

Proof. The convergence follows from ( 13. 5 p since the integrals with the oscillatory functions 
converge to zero. □ 

Below we will need the following lemma: 



Lemma 3.3. Let Conditions SO— S3 hold. Then for any a > 3/2 the bound holds 

supE\\ X (;t)\\h <°° ( 3 - 7 ) 

t>o 

Proof. Denote 

e t (x) := E\x(x,t)\ 2 , x G M 3 . 

The mathematical expectation is finite for almost all x e R 3 by (12.131) with s = —a and the 
Fubini theorem. Moreover, e t (x) = e t for almost all x e R 3 by SI. Formula (13. 5 p implies 

= 7^ /[cos 2 (^)g (A:)-^^(go(A:)P(A:)-P(A;)go(A:)) (3.8) 



(2 



7T 



sin 2 ujt 



P(k)q (k)P(k) 



dk, 



Then et = trg 4 (0) < Ceo- Hence for a > 3/2 we obtain 

E\\x(;t)\\h = e t A 1 + MT^x < C(z/)e 



— (T 



and then (13.71) follows. □ 

We will use also the following result: 

Proposition 3.4. (see J3, Proposition 2.8], J3], Proposition 3.3]). Let Conditions SO— S3 hold. 
Then for any eD, 

Eexp{i(U {t)4j ,<f))} ->exp{-^ Qoo(0,0)}, t oo. (3.9) 

Remark 3.5. In |2J/ the phase space L 2 oc (R 3 ) <g> C 4 has been considered. Nevertheless, all the 
steps of proving the convergence ( 13 .9]) in /#/ remain true if we change L 2 oc (R 3 ) (g> C 4 by L 2 _ a 
with a > 3/2. 



4 Perturbed Dirac equation. 

4.1 Scattering Theory 

To deduce Theorem 12 .91 we construct the dual scattering theory (ll.lOp for finite energy solutions 
using the Boussaid results, [2]. 

Lemma 4.1. (see /H Theorem 1.1]) Let conditions E1-E2 hold and o > 5/2. Then the bound 
holds 

\\P c um\\ L ^ < C(l + |t|)- 3/2 ||^|| L 2, t e R. (4.1) 

Note that for tp e L 2 the solutions Uo(t)if)o and U(t)ipo to problems (12. 12}) and (II. ip . 
respectively, also belong to L? and the charge conservation holds: 

||CWo|| = lhM, \\U (t)i> \\ = (4-2) 
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Here and below || • || is the norm in L 2 . 

For i G 1, introduce the operators U'^it) and U'{t) which are conjugate to the operators 
I7o(*) and U(t) on L 2 : 

%{t)4>) = (U (t)i>, 0), (V, U'{t)4>) = (U(t)ip, 0), t/j, G L 2 . (4.3) 

Here (•, •) stands for the hermitian scalar product in L 2 (IR 3 ,C 4 ). The adjoint groups admit a 
convenient description: 

Lemma 4.2. For G L 2 the function C/o(£)0o = ^>(-, t) is the solution to 

(j){x,t) = [a-V + i/3m](j)(x,t), <p(x, 0) = (f) Q (x). (4.4) 
Proof. Differentiating the first equation of (14. 3 P with ip, G V, we obtain 

The group Uo(t) has the generator 

A = -a ■ V - z/3m. (4.6) 

Therefore, the generator of U' {t) is the conjugate operator 

A' = a • V + i/3m. (4.7) 

Hence, (JO} holds, where 0(t) = ./^(t). 1=1 

Similarly, we obtain 

Lemma 4.3. For G L 2 t/ie function U'(t)<p = 4>{x,t) is the solution to 

<p(x,t) = [a-\7 + i/3m + iV}(f)(x,t), 0(x,O) = 0(x). (4.8) 

Corollary 4.4. i) U' {t) = U (-t), U\t) = U(-t). 
ii) For any G L 2 the uniform bounds hold: 

\Km\ = hi \\u'm\ = Mi *>«• ( 4 - 9 ) 

Under assumptions E1-E2 for U'{t) a bound of type ( 14.1)) also holds: 

\\p c u'm\\ L i„ < c{\ + i*ir 3/2 iHk, t g r (4.io) 

with a > 5/2. 

Now we formulate the scattering theory in the dual representation. 

Theorem 4.5. Let conditions El— E2 and SO— S3 hold and a > 5/2. Then there exist linear 
operators W, r(t) : L 2 — >■ L 2 such that for G L 2 

P c U'{t)(j) = U' Q {t)W(t) + r{t)4), t>0. (4.11) 

and i/je bounds hold 

||r(*)0|| < C{l + ty l / 2 U\\ Li , (4.12) 

£|(Vo,r(t)0>| 2 < C(l + 1)- 1 !!^!!!,, t>0. (4.13) 
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Proof. We apply the Cook method, [HJ Theorem XI. 4]. Fix G L? a and define W<p, formally, 
as 

W<j> = lim U' (-t)P c U'(t)<P = <P+ / U' Q (-r)P c U\r)<pdr. (4.14) 

We have to prove the convergence of the last integral in the norm of L? . First, observe that 

J^(r)0 = ^(r)0, ^U'{r)cj> = AU\r)<j>, r > 
where A' and A' are the generators to the groups Uq(t), U'(t), respectively. Therefore, 

A^(_ T )p c c/'( r )0 = U' g {-t){A' - A^PJJ'ir)^. (4.15) 

We have A' - A' = iV. Furthermore, E2, <KM . (O0j) imply that 

\\%(-t){A-A q )P c U'(t)<P\\ < C \\{A-A g )P c U\t)<P\\=C\\VU\t)<P\\ (4.16) 

< d ||f/'(r)0|| L 2 p < C 2 (l + r)- 3 / 2 ||0|| L 2, r > 0. 

Hence, the convergence of the integral in the right hand side of ( 14. 14ft follows. 
Further, (OH and (041 imply 

r{t)<f> = P c U\t)(j> - U' Q {t)W<j> = -U' Q {t) J —U' Q {-T)P c U'{r)cj>dT. 

Hence (02|) follows by (O) . (Q5|) and (OB) . 

It remains to prove (I4.13p . Applying the Shur lemma we obtain 

E|(^o,r(t)0)| 2 = (g o (x-y),r(t)0(x)®r(t)0(y)) 

< || go || Ll ||r(t)0|| 2 . (4.17) 

Hence, (Q3j) follows by ( ICTj) with p = 1 and ( gJ2) . □ 
4.2 Convergence to equilibrium distribution 

Theorem 12.91 can be derived from Propositions I4.6H4.7I below by using the same arguments as 
in [151 Theorem XII. 5. 2]. 

Proposition 4.6. The family of the measures {P*fit,t G K} is weakly compact in P c Hl^ for 
any e > and o > 5/2 + 5. 

Proposition 4.7. For any G P 

P*iH(<j>) = J exp(z(V,0))P c V t (#) ->exp{-i(2ooW,W^)}, * -> oo. (4.18) 
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Proposition 14.61 provides the existence of the limiting measures of the family P*p,t, and 
Proposition 14.71 provides the uniqueness of the limiting measure, and hence the convergence 
(12. 19(1 . We deduce these propositions with the help of Theorem 14.51 

Proof of Proposition 14.61 First, we prove the bound 

swpE\\P c U{t)tjj \\ n < oo, (4.19) 
t>o 

Representation (I2.14p implies 

P c U(t)ih = P cX (x, t) + P c 0(x, t), (4.20) 
where x( x it) = Uo(t)ip , and <fi(x,t) is the solution to (I2.15p . Therefore, 

E\\P c U{t)MW<E\\P cX mH + E\\P c <P{t)\\ n . (4.21) 

Bound (13 ,7p implies 

snpE\\x(t)\\n<oo. (4.22) 
Further, we have by the Cauchy-Schwartz inequality 

£||(x(0,Ci)C;lk CT < C||OII^_JlOlk^llx(0ll^ < C 3 E\\ X {t)\\ L ^ a = 5/2 + 5 
since the eigenfunctions Q G L 2 S with any s, see Appendix. Therefore 

sup E\\P c x(t)\\n < oo 

t>o 

since P c x(x, t) = x(x, *) - PdXfa t) by (TO]). 

It remains to estimate the second term in the RHS of (14.211) . Choose a Si > such that 
5\ < p — 5 — 5. It is possible due to El. Then the Duhamel representation (12. 151) and bounds 
(gUD and imply 

E\\P e <f>{t)\\n < f E\\P c U{t-s)V X {s)\\ LU/2 _ s ds<C f\l + t - s)- z / 2 E\\Vx{t)\\ L t /2+s ds 

</ J 

< Ci / {\ + t-s)-^E\\ X (S)\\L* ds<C 2 , t>0 (4.23) 

JO 5/2+h-P 

since 5/2 + 5 1 - p < -5/2 - 5. Now (14211 - (OBI imply (091 . 

Now Proposition 14.61 follows from ( 14.191) by Prokhorov theorem [15[ Lemma II. 3.1] as in the 
proof of |15, Theorem XII. 5. 2]. □ 

Proof of Proposition 14.71 We have 



exp(i(^,0))P c ><(#) = J exp(i(P c ip,(j)))p:t(dip) = E expi(P c U(t)ip , 4>) 
Bound (14.131) and Cauchy-Schwartz inequality imply 

\Eexpt(P c U(tW ,<j)) -Eexpi(U (t)i/> ,W<f>)\ = \Eexpi(^ , P c U'{t)4>) - Eexpi(if>o, U' Q {t)W<f>)\ 

< E\(<ifo, r{t)<j>) | < (E(ik, r(t)0) 2 ) 1/2 -> 
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as t — > oo. It remains to prove that 

E^ V i{^U' Q {t)W(j))^e^{-^Q^{W4>,W(j))}, t^oo. (4.24) 

The convergence does not follow directly from Proposition 13.41 since W<p (jL T>. We can approx- 
imate W<p G L 2 by functions from T> since T> is dense in L 2 . Hence, for any e > there exists 
<p e G T> such that 

||W0-0 £ ||<£. (4.25) 

By the triangle inequality 

|£exp i(^ , E^(*W) - exp{-iQ oo (H/0, W^)}| 
< |E exp i<V»o, - £ exp ifa, U' {t)<P £ ) \ 

+E\ expi{U {t)ip ,<f>e) ~ exp{--Q oo (0 £ ,0 e )}| 

+ 1 eM-'QMe, <Pe)} ~ exp{-iQ oo (H/0, W4>)}\. (4.26) 
Let us estimate each term in the RHS of (I4.26p . Theorem 14.51 implies that uniformly in t > 

E\^u'S)(w<p-<p £ ))\ < (E\(^%(t)(w<p- <p £ ))\ 2 ) 1 ' 2 < ||g |li( 2 ||tfo(*)(w^-&)|| 

< C\\W(j)-(j) E \\ < Ce. 

Then the first term is 0(e) uniformly in t > 0. The second term converges to zero as t — > oo 
by Proposition 13.41 since <p £ G T>. At last, the third term is 0(e) by (14. 25ft and the continuity 
of the quadratic form Qoo(0, 4>) in L 2 £g> C 4 . The continuity follows from Corollary 12.81 Now 
convergence (14.241) follows since e > is arbitrary. □ 

5 Appendix: Decay of eigenfunctions 

Here we prove the spatial decay of eigenfunctions. 

Lemma 5.1. Let V satisfy El, and i/j(x) G L 2 (IR 3 ) be an eigenfunction of the Dirac operator 
corresponding to a eigenvalue A G (— m,m), i.e. 

Hip(x) = \ip(x), x G M 3 . 

Then ip G L 2 for all s£l. 

Proof. Denote by Ro(X) = (Hq — A) -1 the resolvent of the free Dirac equation. The equation 
(H + V - X)tp = implies 

^ = R (\)f, where / = —Vip G L\ +p (5.1) 

From the identity 

(-ia ■ V + (3m - X)(ia ■ V - (3m - A) = A - m 2 + A 2 
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it follows that 




(5.2) 



Hence, in the Fourier transform, the first equation of (15.11) reads 



{-a-k + Pm + \)f{k) 
k 2 + m 2 — X 2 



Since |A| < m, we have 




'2+p 



cu\\ H 2 +P < c 1 \\f\\ H ,+ P 




'2+p 



<c 3 U\\ Ll 



Hence, ip G L 2 with any sGMby induction. 



□ 
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